Mon. Not. R. Astron. Soc. OOP. [TlfTTI f 20101 Printed 7 September 2011 (MN WT^ style file v2.2) 



Eurybates — the only asteroid family among Trojans? 



M. Broz^* and J. Rozehnal^'^ 

^Institute of Astronomy, Charles University, Prague, V Holesovickdch 2, 18000 Prague 8, Czech Republic 
^ Stefdnik Observatory, Petfm 205, 11800 Prague, Czech Republic 



Accepted ???. Received ???; in original form ??? 



o 

(N 



o 



> 

O 



0^ 
O 



X 



ABSTRACT 

We study orbital and physical properties of Trojan asteroids of Jupiter. We try to 
discern all families previously discussed in literature, but we conclude there is only 
one significant family among Trojans, namely the cluster around asteroid (3548) Eu- 
rybates. It is the only cluster, which has all of the following characteristics: (i) it is 
clearly concentrated in the proper-element space; (ii) size-frequency distribution is 
different from background asteroids; (iii) we have a reasonable collisional/dynamical 
model of the family. Henceforth, we can consider it as a real collisional family. 

We also report a discovery of a possible family around the asteroid (4709) En- 
nomos, composed mostly of small asteroids. The asteroid (4709) Ennomos is known to 
have a very high albedo pv — 0.15, which may be related to a hypothetical cratering 
event which exposed ice (Fernandez et al. 2003). The relation between the collisional 
family and the exposed surface of the parent body is a unique opportunity to study the 
physics of cratering events. However, more data are needed to confirm the existence 
of this family and its relationship with Ennomos. 
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1 INTRODUCTION 

Trojans of Jupiter, which reside in the neighbourhood of L4 
and I/5 Lagrangian points, serve as an important test of the 
planetary migration theory (Morbidelli et al. 2005). Their 
inclination distribution, namely the large spread of /, can 
be explained as a result of chaotic capture during a brief 
period when Jupiter and Saturn encountered a 1:2 mean- 
motion resonance. Moreover, the Late Heavy Bombardment 
provides the timing of this resonant encounter ~3.8 Gyr ago 
(Gomes et al. 2005). It is thus important to understand the 
population of Trojans accurately. 

There are several unresolved problems regarding Tro- 
jans, however, for example the number of families, which 
is a stringent constraint for collisional models. Roig et al. 
(2008) studied as many as ten suggested families, using rel- 
atively sparse SLOAN data and spectra. They noted most 
families seem to be heterogeneous from the spectroscopic 
point of view, with one exception — the C-type Eurybates 
family. As we argue in this paper, the number of families 
(with parent-body size D > 100 km) is indeed as low as one. 

Another strange fact is the ratio of L4 and I/5 Tro- 
jans. Szabo et al. (2007) used SLOAN colour data to reach 
fainter than orbital catalogues and estimated the ratio to 
N[L4)/N[Lc,) = 1.6 ± 0.1. There is no clear explanation for 
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this, since the chaotic capture as a gravitational interaction 
should be independent of size or L4/I/5 membership. Any 
hypothesis involving collisions would require a relatively re- 
cent disruption of a huge parent body, which is highly un- 
hkely (O'Brien and MorbideUi 2008, D. O'Brien, personal 
communication). This is again related to the actual observed 
number of Trojan families. 

Broz and Vokrouhlicky (2008) studied another resonant 
population, the so called Hilda group in the 3/2 mean- 
motion resonance with Jupiter, and reported only two fami- 
lies: Hilda and Schubart with approximately 200 and 100 km 
parent bodies. This number might be in accord with low col- 
lisional probabilities, assuming the Hilda family is very old 
and experienced the Late Heavy Bombardment (Broz et al. 
2011). 

Levison et al. (2009) compared the observed distribu- 
tion of D-type asteroids and the model of their delivery from 
transneptunian region. They found a good match assuming 
the D-types (presumably of cometary origin) are easy-to- 
disrupt objects (with the strength more than 5 times lower 
than that of solid ice). Note that Trojan asteroids are a mix- 
ture of C- and D-type objects and we have to discriminate 
between them with respect to collisional behaviour. 

All of the works mentioned above are a good motivation 
for us to focus on asteroid families in the Trojan population. 
The paper is organised as follows. First, we describe our 
data sources and methods in Section [2] A detailed study of 
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orbital and physical properties of families (and other 'false' 
groupings) is a matter of Section [3] Section |4] is devoted 
to the modelling of long-term dynamical evolution. Finally, 
there are concluding remarks in Section [5] 



2 METHODS 

2.1 Resonant elements 

We use the symplectic SWIFT integrator (Levison & Dun- 
can 1994) for orbital calculations. Our modifications include 
a second order scheme of Laskar & Robutel (2001) and on- 
line digital filters, which enable us to compute suitable res- 
onant proper elements: libration amplitude d of the a — a' 
oscillations, where a is the osculating semimajor axis of an 
asteroid and a' that of Jupiter, eccentricity e and inclination 
sin/. (In figures, we usually plot a mean value a of semima- 
jor axis plus the libration amplitude d.) We employ their def- 
inition from Milani (1993). The source of initial osculating 
elements is the AstOrb catalogue, version JD = 2455500.5 
(Oct 31st 2010). 

There are actually two independent filters running in 
parallel: in the first one, we sample osculating elements ev- 
ery 1 yr, compute mean elements using the filter sequence 
B, B with decimation factors 3, 3 (refer to Quinn, Tremaine 
& Duncan 1991) a store this data in a buffer spanning 1 kyr. 
We then estimate the libration frequency / by a linear fit of 
(j){t) = A — A' — X, where A, A' are the mean longitudes of an 
asteroid and Jupiter and x = ±60° for L4 or L5 respectively. 
The revolution of angle ^(t) must not be confined to the in- 
terval [0,360°), of course. The amplitude of d is computed 
for the already known / by a discrete Fourier transform. Fi- 
nally, an off-line running-average filter with a window 1 Myr 
is used to smooth the dataQ 

In the second filter, we compute proper eccentricity e 
and proper inclination sin 7 by sampling osculating elements 
(lyr step), computing mean elements using a filter sequence 
A, A, B and decimation factors 10, 10, 3, and then we ap- 
ply a frequency modified Fourier transform (Nesvorny & 
Sidlichovsky 1997), which gives us the relevant proper am- 
plitudes. 

The values of the resonant elements agree very well with 
those listed in the AstDyS catalogue by Knezevic & Milani 
(see Figure [T|. There are only few outliers, probably due to 
a different time span of integration. We computed proper 
elements for 2647 1/4 and 1496 I/5 Trojan asteroids|fj This 
sample is roughly twice larger than previously analysed. The 
ratio of populations valid for H <Vb mag asteroids is thus 
N{U)/N{L^) ~ 1.8. 

The overall distribution of Trojans in the (d, e,sin/) 
space is shown in Figure [21 Note there is only one clus- 
ter visible immediately in the bottom-left panel — around 
(3548) Eurybates. The reason is its tight confinement in in- 
clinations (sin / = 0.125 to 0.135). 



^ Equivalently, we may compute the amplitude D of mean lon- 
gitudes A — A'. Anyway, there is a linear relation between d and 
D. 

^ The data are available in an electronic form on our web site 
|http : //sirrah. troja.mff ■ cuni . cz/~mlra/mp/ [ We use also one- 
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Figure 1. Comparison of the resonant eccentricity calculated by 
our code to that of Knezevic & Milani (AstDyS catalogue). There 
is a line x = y to aid a comparison. 

2.2 Hierarchical clustering 

In order to detect clusters in the resonant element space we 
use a hierarchical clustering method (Zappala et al. 1994) 
with a standard metric d\ , with 5a substituted by d. We run 
the HCM code many times with various starting bodies and 
different cut-off velocities Ucutoff and determine the num- 
ber of bodies A'^ in the given cluster. We find the A'^(wcutoff) 
dependence a very useful diagnostic tool. We can see these 
dependences for L4 and is Trojans in Figure [3] 

It is easy to recognise, if a cluster has a concentration 
towards the centre — even at low Wcutoff it must have more 
than one member (A'^ 3> 1). It is also instructive to compare 
clusters with a random background (thin lines), which we 
generated artificially by a random-number generator in the 
same volume of the (d, e, /) space. Insignificant (random) 
clusters usually exhibit an abrupt increase of Af at a high 
cut-off velocity. 

As starting bodies we selected those listed in Roig et al. 
(2008). Only three clusters, namely the Eurybates, Aneas, 
1988 RGio seem to be somewhat concentrated, i.e., denser 
than the background. The Hektor cluster is also concen- 
trated but it contains only a relatively low number of mem- 
bers (20 to 70) before it merges with the background. In 
other words, smaller asteroids do not seem concentrated 
around (624) Hektor. Remaining clusters are more or less 
comparable to the background. 

Nevertheless, we report a detection of a previously un- 
known cluster around (4709) Ennomos in L5 . It is relatively 
compact, since the minimum cut-off velocity is 70m/s only. 
The cluster contains mostly small bodies which were discov- 
ered only recently. 

Finally, let us point out a very tight cluster around 
(9799) 1996 RJ, associated already at Wcutoft = 20m/s. It 
is located at high inclinations and contains 9 bodies, three 
of them having short arcs. The cluster seems peculiar in the 
osculating element space too since it exhibits a non-random 
distribution of nodes and perihelia (see Table [1} . This is 
similar to very young families like the Datura (Nesvorny et 



apparition orbits for the purposes of physical studies. Of course, 
orbital studies require more precise multi-apparition data. 
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Figure 2. The resonant elements (a = a + d, sin/) and (e,sin/) for L4 and L5 Trojans. The crosses indicate relative sizes of bodies, 
taken either from the AstOrb catalogue or computed from absolute magnitude H and geometric albedo py. In this plot, we assumed 
pv = 0.058 for L4 Trojans and 0.045 for those in L5 (it corresponds to medians of known pv's). The asteroids (3548) Eurybates in L4 
and (4709) Ennomos in L5, around which significant clusters are visible, are shown in red. Moreover, the asteroid (9799) 1996 RJ in L4, 
which is surrounded by a small cluster, is denoted by a blue circle. (This cluster is so tight, that its members are located inside the circle 
on the (e, sin/) plot.) 



al. 2006) and it makes the 1996 RJ cluster a particularly in- 
teresting case with respect to coUisional dynamics. Because 
one has to use slightly different methods for studies of such 
young families we postpone a detailed analysis to a next 
paper. 

Let us compare Trojan clusters to the well known aster- 
oid families in the outer Main Belt (Figure [4]). Most families 
(e.g., Themis, Koronis, Eos) exhibit a steady increase of 
until they merge with another families or the entire outer 
Main Belt. Eurybates, Aneas and 1988 RGio are the only 
Trojan clusters which behave in a similar fashion. The Ver- 
itas family (dynamically young, Nesvorny et al. 2003) ex- 
hibits a different behaviour — for a large interval of Wcutofi 
the number of members A'^ remains almost the same, which 
indicates a clear separation from the background population. 
With respect to the A''(iicutofi) dependence, the Ennomos 
cluster is similar to Veritas. 

2.3 Size-frequency distribution 

At first, let us assume a single value of albedo for all family 
members. It is a reasonable assumption if the family is of 
collisional origin. We can then calculate sizes from absolute 
magnitudes and construct size-frequency distributions. Fig- 
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Figure 4. The A'^('Ucutofr) dependence for seven outer main-belt 
families. If we would consider only a subset of asteroids brighter 
than H = 15 mag, which is an approximate observational limit 
for Trojans, the N{vcutoa) dependencies would be qualitatively 
the same, only slightly shifted to larger cut-off velocities. 

ure[5]shows a comparison of SFD's for the clusters detected 
by the HCA/jf] and for the whole population of L4 and L5 
Trojans. 

^ Of course, we have to select a 'suitable' value of the cut-off 
velocity for all clusters. Usually, we select that value for which 
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Table 1. List of nine members of the (9799) 1996 RJ cluster and their proper (a, e,sin/) and osculating (f!osc,roosc) elements and 
absolute magnitude H. Note the distribution of nodes and perihelia is not entirely uniform. Asteroids with short-arc orbits (<60 days) 
are denoted by a * symbol. 



number 


designation 


a 


e 


sin / 






H/mag 


9799 


1996 RJ 


5.2252 


0.0412 


0.5269 


115.4 


259.6 


9.9 


89938 


2002 FR4 


5.2324 


0.0394 


0.5274 


70.0 


23.1 


12.5 


226027 


2002 EKi27 


5.2316 


0.0399 


0.5263 


62.8 


352.9 


12.6 


243316 


2008 RL32 


5.2340 


0.0398 


0.5268 


27.3 


358.2 


12.8 




2005 MG24 


5.2275 


0.0404 


0.5252 


172.3 


236.5 


13.1 




2008 OW22 


* 5.2276 


0.0401 


0.5274 


53.7 


340.9 


13.9 




2009 RA17 


* 5.2258 


0.0409 


0.5272 


257.7 


194.5 


13.7 




2009 RK63 


* 5.2305 


0.0407 


0.5260 


56.4 


5.6 


12.8 




2009 SR30 


5.2362 


0.0409 


0.5258 


103.6 


22.0 


13.3 
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Figure 3. Left panel: The dependence of the number of family members N on the cut— off velocity t'cutotr computed by the hierarchical 
clustering method. Only clusters among L4 Trojans are included in this plot. Middle panel: The same A'^(f cutoff) dependence for L5 
Trojans. Right panel: Artificial clusters selected from random distribution of asteroids generated in the same volume of the (d, e, sin/) 
space. 



A slope 7 of the cumulative distribution N{>D) oc 
is an indicative parameter. For L4 and L5 Trojans it equals 
to — 2.0 ± 0.1 and — 1.9 ± 0.1 in the intermediate size range 
15 to 60 km. (These numbers are compatible with the study 
of Yoshida & Nakamura 2008.) The slope is steeper at large 
sizes. The uncertainties are mainly due to a freedom in se- 
lection of the size range and the difference between L4 and 
L5 SFD's does not seem significant. The clusters have typ- 
ically similar slope as background (within 0.1 uncertainty), 
thought sometimes the results are inconclusive due to small 
number of members. The Eurybates family with —2.5 ±0.1 
slope is on the other hand significantly steeper than the 
mean slope of the whole Trojan population^ There are two 
more groups which exhibit a relatively steep slope, namely 
Laertes in L4 (7 = -3.1) and 1988 RGio in L5 (7 = -2.6). 

We should be aware, however, that even the background 
exhibits a trend with respect to inclinations (see Figure 
Slope 7 typically decreases with inclination sin/, which is 
especially prominent in case of the L4 cloud. We have to ad- 
mit if we compare the Eurybates family to its surroundings 



^("cutoff) is fiat. Size-frequency distribution is not very sensitive 
to this selection anyway. 

* Thought the number of the Eurybates members (105) is so 
small that it almost does not affect the mean slope of the whole 
L4 population. 



only (sin / = 0.1 to 0.15), the difference in slopes is not so 
prominent. An interesting feature of the L5 cloud is a dip 
in the interval sin / = 0.05 to 0.1. This corresponds to the 
approximate location of the 1988 RGio group. 

The 7(sin/) dependence among Trojans is not unique. 
E.g. low-inclination bodies in the J3/2 resonance also have 
the SFD steeper than background (7 = —2.5 ± 0.1 versus 
— 1.7 ±0.1), without any clear family and a few big interlop- 
ers. May be, this feature reflects different source reservoirs of 
low- and high-inclination bodies among Trojans and 
It may be also in concert with a colour-inclination depen- 
dence reported by Szabo et al. (2007). 

We also test albedo distributions dependent on size, 
since the measurements by Fernandez et al. (2009) sug- 
gested small Trojans are signiflcantly brighter and thus 
smaller. Large asteroids hayepv = 0.044 ±0.008 while small 
pv = 0.12 ± 0.06. This is a significant change of the SFD, 
which occurs around the size D ~ 30 km. The SFD thus be- 
comes shallower below this size, e.g. for Eurybates we would 
have 7 = —1.6 and for L4 Trojans 7 = —1.5, so the SFD's 
become comparable with respect to the slope. Thought, as 
we stated above, for a real coUisional family we expect the 

^ Both Trojan and J3/2 regions are dynamically unstable during 
Jupiter— Saturn 1:2 mean motion resonance, so we expect the same 
bodies entering Trojan region may enter J3/2. 
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Figure 5. Left panel: size distributions of L4 Trojans and the following clusters (there is a selected cut— off velocity in the parenthesis): 
Eurybates (^cutoff = 50m/s), Laertes (94), Hektor (160), Teucer (175), Sinon (163), 1986 WD (120). Right panel: SFD's of L5 Trojans 
and the following clusters: 1988 RGio (at i^cutofl = 130 m/s), Aneas (150), Asios (155), Panthoos (130), Polydoros (130). 
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Figure 6. Slopes 7 of the size-frequency distributions N{>D) 
for L4 and L5 Trojans and their dependence on the inclina- 
tion sin/. The range of diameters for which the SFD's were fitted 
is Cmin = 12 km, Dmax = 30 km. Thin lines were calculated 
for different ranges, which were varied as D^in G (10, 15) km, 
fmax G (20, 40) km. Their spread indicates the uncertainty of 7 
in a given interval of sin/. The populations are observationally 
complete down to /? ~ 10 km, because the characteristic change 
of slope due to incompleteness occurs at smaller sizes (see also 
Yoshida and Nakamura 2008). 



albedo distribution to be rather homogeneous and indepen- 
dent of size. 



2.4 Colour and spectral data 

We used the Sloan Digital Sky Survey Moving Object cat- 
alogue version 4 (SDSS-M0C4) to check the families are 
spectrally homogeneous, as we expect. Due to a larger un- 
certainty in the u colour in SDSS-MOC4, we used the color 



indices a* and i — z, where a* = 0.89((; — r)-f 0.45(r — i) — 0.57 
(defined by Parker et al. 2008). 

The result is shown in Figure [T] It is clearly visible 
that the distribution of the Eurybates family in the space of 
(a*,i — z) colours is different from the Trojan background. 
On contrary, the 1988 RGio group covers essentially the 
same area as the background. The Aneas is only slightly 
shifted towards larger a* and i~ z with respect to the back- 
ground. There is a lack of data for the Ennomos group — 
three bodies are not sufficient to compare the colour distri- 
butions. 

Alternatively, we may use principal component analysis 
of the SDSS colour indices. We use only data with uncertain- 
ties smaller than 0.2 mag, which resulted in 70 887 records. 
We calculated eigenvalues (Ai, 2,3,4 = 0.173, 0.0532, 0.0249, 
0.0095), corresponding eigenvectors and constructed the fol- 
lowing three principal components (Trojanova 2010): 

PCi = 0.235 (ti - ff) + 0.416 (g-r) + 0.598 (3 -i) 

+ 0.643 (5 -^), (1) 

PC2 = 0.968 (m -3)- 0.173 (g-r-)- 0.147 (s-i) 

- 0.106 (5 -z), (2) 

PCs = 0.078 {u-g) + 0.601 {g - r) + 0.330 {g - i) 

-0.724(5-2), (3) 

which have a clear physical interpretation: PCi corresponds 
to an overall slope, PC2 is a variability in the u band, and 
PC3 a depth of the 1 /im absorption band. The Eurybates 
family is different from Trojans in all three principal compo- 
nents (mean PCi of the Eurybates members is smaller, PC2 
and PCs larger). The Aneas group has the same distribu- 
tion of PC2 and PCs as Trojans and the 1988 RGio group 
is similar to Trojans even in all three components. 

Hence, we confirm the Eurybates family seems distinct 
in color even in the fourth version of the SDSS-MOC. This 
fact is consistent with previous work of Roig et al. (2008), 
who used third version of the same catalogue and classified 
Eurybates family members as C-type asteroids. 

Finally, note that De Luise et al. (2010) pointed out an 
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Figure 7. Left panel: The {a*,i — z) colours for the L4 Trojans 
(gray dots) and the Eurybates family (black dots with error bars). 
The distributions differ significantly in this case. Right panel: A 
similar comparison for the L5 Trojans and the 1988 RGio group, 
which seem to be indistinguishable. 



absence of ice spectral features at 1.5 and 2.0 ^im on several 
Eurybates members and Yang and Jewitt (2007) concluded 
the same for (4709) Ennomos. This puzzling fact may in- 
dicate that pure ice covers at most 10 % of the Ennomos 
surface. 



2.5 Impact disruption model 

We use a simple model of an isotropic disruption from 
Farinella et al. (1994). The distribution of velocities "at in- 
finity" follows the function 



dN{v) = Cv{v^ +vL) 



(4) 



with the exponent a being a free parameter, C a normali- 
sation constant and «osc the escape velocity from the parent 
body, which is determined by its size Rpb and mean den- 
sity ppB- The distribution is cut at a selected maximum 
allowed velocity iimax to prevent outliers. We typically use 
'^max = 300 m/s. The orientations of velocity vectors in space 
are assigned randomly. We assume the velocity of fragments 
is independent on their sizelfl 

There are several more free parameters, which deter- 
mine the initial shape of the family in the space of proper 
elements: initial osculating eccentricity ei of the parent body, 
initial inclination ii, as well as true anomaly /imp and argu- 
ment of perihelion ujimp at the time of impact disruption. 

An example of a synthetic family just after disrup- 
tion and its comparison to the observed Eurybates family 
is shown in Figure [S] Usually, there is a significant disagree- 
ment between this simple model of impact disruption and 
the observations. Synthetic families usually look like thin 
'filaments' in the (d, e, sin J) space, which are curved due to 
the mapping from osculating elements to resonant ones. On 
the other hand, observed groups among Trojans are much 
more spread. However, this only indicates an importance of 



^ If we use a size-dependent relation for velocities similar to 
Vokrouhlicky et al. (2006), our results do not change much, be- 
cause the overall shape of the velocity distribution is quite similar 
to the size-independent case. 



further long-term evolution by chaotic diffusion and possibly 
by planetary migration^ 

In case of the Ennomos group members, they are dis- 
tributed mostly on larger semimajor axes than (4709) En- 
nomos, thought isotropic impact disruptions produce frag- 
ments distributed evenly on larger and smaller a. May be, 
it is an indication of an anisotropic velocity field? Or a dif- 
ferent parent body of this cluster? 



2.6 Planetary migration 

If asteroid families are very old, planetary migration might 
influence their current shape. In order to study of late stages 
of planetary migration, which is caused by interactions with 
a planetesimal disk, we construct the following model. We 
treat the migration analytically within a modified version of 
the numerical symplectic SWIFT-RMVS3 integrator (Levi- 
son & Duncan 1994), which accounts for gravitational per- 
turbations of the Sun and four giant planets and includes 
also an energy-dissipation term, as described in Broz et al. 
(2011). The speed of migration is characterised by the ex- 
ponential time scale Tmig and the required total change of 
semimajor axis Oi — Sf . We use an eccentricity damping for- 
mula too, which simulates the effects of dynamical friction 
and prevent an unrealistic increase of eccentricities (Mor- 
bidelli et al. 2010). The amount of damping is determined 
by the parameter edamp. 

We try to adjust initial orbital parameters of planets 
and the parameters of migration in such a way to end up 
at currently observed orbits. The integration time step is 
At = 36.525 days and the time span is usually equal to 
3rniig, when planetary orbits almost stop to migrate. 



2.7 Inefficient Yarkovsky/YORP effect 

On long time scales, the Yarkovsky thermal force might 
cause significant perturbations of orbits. We use an imple- 
mentation of the Yarkovsky thermal effect in the SWIFT 
N-body integrator (Broz 2006). It includes both the diurnal 
and seasonal variants. 

The YORP effect (thermal torques affecting spin states; 
Vokrouhlicky et al. 2006) was not taken into account in 
our simulations. The reason is that the respective time 
scale tyorp is of the order 100 My to 1 Gyr So as a 'zero' 
approximation, we neglect the YORP effect on these "short" 
time scales and keep the orientations of the spin axes fixed. 

For Trojan asteroids captured in a zero-order mean mo- 
tion resonance, the Yarkovsky perturbation only affects the 
position of libration centre (Moldovan et al. 2010). Note 
that the perturbation acts 'instantly' — there is no system- 
atic secular drift in eccentricity nor in other proper elements 
which is an important difference from first-order resonances, 
where a e-drift is expected (Broz & Vokrouhlicky 2008, Ap- 
pendix A). This is another reason we do not need a detailed 
YORP model here. 

The thermal parameter we use are reasonable estimates 
for C/X-type bodies: psurf = Pbuik = 1300 kg/m^ for the 
surface and bulk densities, K = 0.01 W/m/K for the surface 



^ Only very young clusters like the Karin family (Nesvorny et al. 
2002) exhibit this kind of a 'filament' shape. 
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Figure 8. A comparison between the observed Eurybates family (open circles) and synthetic families (crosses) just after the impact 
disruption computed for several values of /imp = 0° , 45° , 90° , 135° , 180° and wimp = 30° , Rpb = 47 km, ppb = 1300 kg/m^ . Different 
geometry in /, uj produces a slightly different cluster, nevertheless, it is always tighter than the observed family. The position of the 
asteroid (3548) Eurybates is denoted by a square. 



thermal conductivity, C = 680J/kg for the heat capacity, 
A = 0.02 for the Bond albedo and eiR = 0.95 for the thermal 
emissivity. 



3 ASTEROID FAMILIES AND 
INSIGNIFICANT GROUPS 

In this section, we briefly discuss properties of selected clus- 
ters: Eurybates, Ennomos, Aneas and 1988 RGio. We focus 
on these four clusters, since they seem most prominent ac- 
cording to our previous analysis. 



3.1 Eurybates family 

The Eurybates family can be detected by the hierar- 
chical clustering method for cut-off velocities i^cutoff ~ 
38 to 78m/s, when it merges with Menelaus (see Figure |3]). 
Yet, we do not rely just on the HCM! Another selection 
criterion we use is a 'meaningful' shape of the family and 
its changes with respect to Wcutoff- A very important char- 
acteristic of the Eurybates family at low values of Ucutoff 
is a tight confinement of inclinations (sin / within 0.01). It 
breaks down at «cutoff — 68m/s, so we consider this value 
as an upper limit. The Eurybates family is also confined 
in semimajor axis, being approximately twice smaller than 
other groups. 

The diameter of the parent body is DpB = 97 km for 
albedo pv = 0.055 if we sum the volumes of the known 
bodies. Of course, in reality it is slightly larger due to obser- 
vational incompleteness. If we prolong the slope of the SFD 
7 — —2.5 down to zero we obtain DpB = 110 km. The ge- 
ometric method of Tanga et al. (1999) gives an upper limit 
DpB ~ 130 km. 

Spectral slopes of family members are rather homoge- 
neous and correspond to C/P- types (Roig et al. 2008). 



3.2 Ennomos group 

The cluster around (4709) Ennomos can be recognised for 
a wide interval of cut-off velocities Wcutoff G (69, 129) m/s 
when it stays compact and confined in inclinations (sin / = 
0.451 to 0.466). Very probably, there are several interlopers, 
because we can count 4 to 10 asteroids in the surroundings. 



i.e., in the same volume of the (d, e,sin/) space (see Fig- 
ure [9|. Since small bodies dominate the Ennomos group we 
suspect large bodies might be actually interlopers. 

A very intriguing feature is a high albedo of (4709) En- 
nomos pv — 0.15 measured by Fernandez et al. (2003). 
Apart from other explanations, the authors speculated it 
may result from a recent impact which covered the surface 
with pristine ice. If true the relation between the fresh sur- 
face and the coUisional family might be a unique opportunity 
to study cratering events. 

We cannot exclude a possibility that (4709) Ennomos 
is actually an interloper and the family is not related to 
it at all. Nevertheless, our hypothesis is testable; family 
members should exhibit a similarity in spectra and albe- 
dos. The only information we have to date are SDSS colours 
for three members: 98362, 2005 YG204 are probably C-types 
and 2005 AR72 is a D-type. In case new data become avail- 
able we can remove interlopers from our sample and improve 
our analysis. 

The size distribution of the Ennomos group is barely 
constrained, since small bodies are at the observational 
limit. Moreover, removal of interlopers can change the SFD 
slope completely (from 7 = —1.4 to —3.2 or so). The mini- 
mum parent body size is about Dpp — 67 km if all members 
have high albedo pv ~ 0.15. 

3.3 Group denoted Aneas 

The Aneas group looks like a middle portion of the L5 cloud 
with approximately background density. It spans whole 
range of semimajor axes, as background asteroids do. 

The minimum size of a hypothetical parent body is 
DpB = 160 to 170km (for albedo pv = 0.055 to 0.041). 
This size is very large and an impact disruption of such 
body is less probable (see Section [4.4|) . The size-frequency 
distribution is shallow, with approximately the same slope 
as background. 

According to Roig et al. (2008) the colours are rather 
homogeneous and correspond to D-types, with ~ 10 % of 
probable interlopers. 

3.4 Group denoted 1988 RGio 

The group around asteroid (11487) 1988 RGio again looks 
like a lower portion of the L5 cloud at low inclinations, with 
sin / G (0.06,0.1). The SFD is steeper (7 = -2.6 ±0.1) than 
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Figure 9. A detail of the L5 Trojan population where the Ennomos group is visible. Left panel: resonant semimajor axis a vs eccentricity e. 
Only asteroids occupying the same range of inclinations as the Ennomos group sin / £ (0.448, 0.468) are plotted to facilitate a comparison 
with the density of surroundings space (background). The sizes of plus signs are proportional to diameters of the asteroids. Probable 
family members are denoted by small red circles and probable interlopers by small grey crosses. Right panel: a vs inclination sin I, with 
range of eccentricities e S (0.02, 0.045). 



surroundings in L5 and the resulting parent body size D ~ 
60 km is relatively small. The colours seems heterogeneous 
(Roig et al. 2008) and we can confirm this statement based 
on the new SDSS-MOC version 4 data. 

The remaining clusters (Hektor, Teucer, Sinon, 
1986 WD, Laertes, Asios, Polydoros, Panthoos, etc.) may be 
characterised as follows: (i) they have a density in (d, e, sin J) 
space comparable to that of background (surroundings); 
(ii) when identified by the HCM their semimajor axes span 
the whole range of Trojan region; (iii) the slopes of their 
SFD's are comparable to the background; (iv) they are often 
inhomogeneous with respect to colours (according to Roig 
et al. 2008). These reasons lead us to a conclusion that these 
clusters are not necessarily real coUisional families. 



4 LONG-TERM EVOLUTION OF TROJAN 
FAMILIES 

4.1 Evolution due to chaotic diffusion 

We try to model long-term evolution of the Eurybates fam- 
ily. At first, we generate a synthetic family (consisting of 42 
bodies) by an impact disruption of the parent body with 
required size. Then we integrate the synthetic family and 
compare it at particular time to the observed Eurybates 
family. The time span of the integration is 4Gyr. 

The main driving mechanism is slow chaotic diffusion 
(the Yarkovsky effect is present but inefficient in the Trojan 
region). Initially, the spread of inclinations of the synthetic 
family is consistent with the observed one. On the other 
hand, the shape in (a, e) elements is clearly inconsistent. 

Since the inclinations evolve only barely we focus on the 
evolution of in the (a, e) plane (see Figure [10)) . The point is 
the synthetic family, namely the 'filament' structure, has to 
disperse sufficiently. After 500 Myr it is still recognisable but 
after 1 Gyr of evolution it is not. So we may constrain the 
age of the Eurybates family from 1 to 4 Gyr|j 



A similar analysis for the Ennomos group indicates that 
chaotic diffusion is faster in this region (given the large in- 
clination) and the most probable age thus seems to be from 
1 to 2 Gyr. Beyond 2 Gyr the inclinations of the synthetic 
family become too large compared to the observed Ennomos 
group, while the eccentricites are still compatible. 

We try to model Aneas and 1988 RGio groups too (see 
FigurefTTjl. In these two cases, there is a strong disagreement 
between our model and observations. The observed groups 
are much larger and chaotic diffusion in respective regions is 
very slow. Even after 4 Gyr of orbital evolution, the synthetic 
family remains too small. 

The only free parameter which may substantially 
change our results is the initial velocity distribution. 
Theoretically, the distribution might have been strongly 
anisotropic. However, we cannot choose initial velocities en- 
tirely freely, since their magnitude should be comparable to 
the escape velocity from the parent body, which is fixed by 
the size Dps and only weakly dependent on a-priori un- 
known density ppB- 

Another solution of this problem is possible if we assume 
families are very old and they experienced perturbations due 
to planetary migration. 

4.2 Stability during planetary migration 

The major perturbation acting on Trojans are secondary res- 
onances between the libration period Pji/i of the asteroid 
in the J 1/1 mean- motion resonance with Jupiter and the pe- 
riod P1J-2S of the critical argument of Jupiter-Saturn 1:2 
resonance (Morbidelli et al. 2005) 

-Pji/i = nPij_2S , (5) 

where n is a small integer number. Typical libration periods 
are Pji/i ~ 150 yr and Pij_2S changes as planets migrate 
(it decreases because Jupiter and Saturn recede from their 
mutual 1:2 resonance) 



^ We verified these estimates by a 2-dimensional Kolmogorov- 
Smirnov test of the (a, e) distributions: initially the KS distance 
is Oks = 0.30 and the probability PKs(>^) = 0.02, which means 
the distribution are incompatible. At t = 1 Gyr, the values are 



^KS = 0.20 and Pks(>^) = 0.32, which indicates a reasonable 
match. 

^ Another source of instability might be a secondary resonance 
with P2.J— 5S (the so called Great Inequality period) thought it is 
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Figure 10. Orbital evolution of the synthetic family and its comparison with the observed Eurybates family. Left panel: the situation 
in the (a, e) plane at 500 Myr. Middle panel; the situation after 4 Gyr. Chaotic diffusion disperses the synthetic family in course of time 
(see shaded tracks of particles). Right panel: the (a, sin/) plane at the same time. Inclinations evolve only barely. 
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Figure 11. Evolution of the synthetic family over 4 Gyr versus 
the observed Aneas group. Chaotic diffusion is slow and it seems 
impossible to match the large spread of the observed group even 
after 4 Gyr. 
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All synthetic families are strongly unstable when 
Pij_2s — 150 yr and even during later stages of migration 
with Pij_2S — 75 yr the eccentricities of family members are 
perturbed too much to match the observed families like Eu- 
rybates or Ennomos (see Figure[T2}. There are practically no 
plausible migration scenarios - regardless of time scale Tmig 
- which would produce a sufficiently compact group, unless 
Jupiter and Saturn are almost on their current orbits. We 
tested Tmig ~ 0.3, 3, 30 Myr and even for Aaj = ajf — aji as 
small as —0.08 AU and Aas — +0.25 AU the perturbation 
was too strong. The reason is that one has to avoid n = 2 
secondary resonance to preserve a low spread of a synthetic 
family. 

Let us conclude if any of Trojan families was created 
during planetary migration and if the migration was smooth 
(exponential) the family cannot be visible today. However, 
we cannot exclude a possibility that final stages of migration 
were entirely different, e.g., similar to the 'jumping- Jupiter' 
scenario (Morbidelli et al. 2010). 



weaker than Pij_2S- We find no asteroids perturbed by secondary 
resonances connected with P^j-rs or P4j_9s which are present 
'en route'. Neither Uranus nor Neptune play an important role. 



t [Myr] 

Figure 12. Evolution of a synthetic family during late phases of 
planetary migration (rmig = 30 Myr in this case). Top panel: the 
state at Myr, middle: 100 Myr, bottom: the respective orbital 
evolution of Jupiter and Saturn. The family is almost destroyed 
and it is definitely incompatible with the observed Eurybates fam- 

iiy- 



4.3 Families lost by the ejection of fragment 
outside the resonance 

We studied a possibility that some families cannot be identi- 
fied because the breakup occurred at the outskirts of the sta- 
ble libration zone and some fragments were ejected outside 
the J 1/1 resonance. We thus chose 30 largest asteroids near 
the edge of the L4, libration zone and we simulated breakups 
of these asteroids which create families with 30 fragments 
each. We assumed the diameter of all parent bodies equal 
to DpB ~ 100km and their density ppB = 1.3 gem"''. The 
breakups always occurred at the same geometry /imp = O'', 
i^imp = 30°. After the breakup, we calculated proper ele- 
ments of the family members and plotted their distribution 
(see Figure [13} . We can see all 30 synthetic families can be 
easily identified. In most cases, more than 95 % of family 
members remained within the stable libration zone. We can 
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Figure 13. Proper eccentricities and inclinations of 30 synthetic 
families (black dots), which originated near the border of stable 
libration zone, compared to the observed L4 Trojans (gray dots). 

thus conclude that the ejection of fragments outside the li- 
bration zone does not affect the number of observed families 
among Trojans. 

4.4 Collisional rates 

We can estimate collisional activity by means of a simple sta- 
tionary model. Trojan-Trojan collisions play a major role 
here, because Trojans are detached from the Main Belt. 
In case of Eurybates, the target (parent body) diameter 
ftarget = 110 km, the mean impact velocity Mmp = 4.7km/s 
(DeirOro et al. 1998), the strength QJ, = 10^ J/kg (Benz & 
Asphaug 1999) and thus the necessary impactor size (Bottke 
et al. 2005) 

rfdisrupt = (SQ^/V^mp)^''^ Aargot ~ 23 km . (6) 

Number of ^23 km projectiles among L4 Trojans is 
?^projcct ~ 371 and we have ntargct ~ 8 available targets. 
An intrinsic collision probability for Trojan-Trojan colli- 
sions Pi = 7.8 X 10"^*km"Vr"^ (Dell'Oro et al. 1998) and 
corresponding frequency of disruptions is 

r rj -^target „ in— H —1 /'yX 

/disrupt — J i Jlprojcct Wtargot ~ ( ■ iU yt . (I) 

Over the age of the Solar System Tss — 4 Gyr (after the 
LHB), we have a very low number of such events nevents = 
Tss/disrupt — 0.28. This number seems to be in concert with 
only one D 100 km family currently observed among Tro- 
jans0 In a less-likely case, the material of the Eurybates 
parent body was very weak and its strength may be at most 
one order of magnitude lower, QJj ~ 10* J/kg (see Lein- 
hardt & Stewart 2009, Bottke et al. 2010). We then obtain 
^disrupt — 10 km and rieventa — 1.0, SO the conclusion about 
the low number of expected Trojan families remains essen- 
tially the same. 

The parent body of Aneas group is 1.5 larger and 
consequently the resulting number of events is more than 

A similar stationary estimate valid for the Main Asteroid Belt 
gives the number of events 12 while the number of observed fam- 
ilies with DpB > 100 km is about 20 (Durda et al. 2007). These 
two numbers are comparable at least to order-of-magnitude. 



one order of magnitude lower. On the other hand, clusters 
with smaller parent bodies (Dpb <^ 100 km) or significantly 
weaker (Q^ <C 10^ J/kg) might be more frequent. 

During the Late Heavy Bombardment epoch we may as- 
sume a substantial increase of collisional activity (Levison et 
al. 2009). Hypothetical old families were however probably 
'erased' due to the late phases of planetary migration (see 
Section [4.21) unless the migration time scale for Jupiter and 
Saturn was significantly shorter than the time scale of the 
impactor flux from transneptunian region which is mainly 
controlled by the migration of Uranus and Neptune. 



5 CONCLUSIONS 

Increasing number of Trojan asteroids with available proper 
elements enables us to get new insights into this impor- 
tant population. Essentially, new faint /small asteroids filled 
the 'gaps' in the proper-element space between previously 
known clusters and today it seems most clusters are rather 
comparable to background. One should be aware that the 
number of families among Trojans may be low and one 
should not take the number of ~ 10 families as a rule. 

Only the C-type Eurybates family fulfils all criteria to 
be considered a collisional family. This is probably also true 
for the newly discovered Ennomos group. Moreover, there 
might be a potentially interesting relation between the high- 
albedo surface of (4709) Ennomos and the collisional family 
thought we do not have enough data yet to prove it inde- 
pendently (by colours, spectra or albedos). 

Note there may exist clusters among Trojans which are 
not of collisional origin. They may be caused by: (i) dif- 
ferences in chaotic diffusion rates; (ii) a/e//-dependent effi- 
ciency of original capture mechanism; or (iii) it may some- 
how reflect orbital distribution in source regions. 

We cannot exclude a hypothetical existence of old fam- 
ilies which were totally dispersed by dynamical processes, 
e.g., by perturbations due to planetary migration which is 
especially efficient in the Trojan region. 

Finally, note there seem to be no D-type families any- 
where in the Solar System — neither in the Trojan re- 
gion, nor in the J3/2 (Broz et al. 2011) and Cybele regions 
(Vokrouhlicky et al. 2010). May be, the D-type parent bod- 
ies are too weak and the target is completely pulverized 
during a collision? This might have important implications 
for collisional models of icy bodies. 
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